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ff~j | We consider pair production and decay of fundamental unstable particles 

in the framework of a modified perturbation theory (MPT) treating resonant 
contributions in the sense of distributions. The cross-section of the process is 
O-f calculated within the NNLO of the MPT in a model that admits exact solu- 

tion. Universal massless-particles contributions are taken into consideration. 
A comparison of the outcomes with the exact solution demonstrates excellent 
convergence of the MPT series at the energies near and above the maximum of 
the cross-section. 
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A description of the processes of productions and decays of fundamental unstable 
particles to satisfy the up-to-date requirements must provide, on the one hand, gauge 
cancellations and unitarity and, on the other hand, enough high accuracy of calcula- 
tion of resonant contributions of unstable particles. Unfortunately, in the framework 
of conventional perturbation theory (PT) a simultaneous fulfilling of these require- 
ments is obstructed by divergences caused by resonant contributions. For this reason 
in the propagators of unstable particles the Dyson resummation is usually applied, 
which shifts the resonant singularities out of the region of physical momenta. How- 
ever, a resummation mixes the PT orders, which generally leads to violation of the 
gauge cancellations. So simultaneously with using the Dyson resummation an appli- 
cation of additional tricks is required. 

Among various approaches that include such tricks, the most known one is based 
on the Laurent expansion of the amplitude around the complex poles of the resonant 
propagators. Each term of this expansion is considered expanded in the framework of 
the conventional PT, as well, but a certain portion of the self-energy is not involved 
in the latter expansion as having been absorbed by the shift of the point of singularity 
(the remnant of the Dyson resummation). The gauge cancellations are completely 
maintained in this approach. However, the precision of the description vastly falls 
at the increasing of a distance from the resonant region, and an uncertainty arises 
at calculating the residues in the complex-poles. Nevertheless, in the vicinity of the 
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resonant region the pole expansion in many cases is suitable for applications. In par- 
ticular, at LEP2 the loop corrections to the W^-pair production were calculated in 
the double pole approximation (DPA) [TJ [2] , the leading approximation in the pole 
expansion. Unfortunately, at international linear collider (ILC) [3] the accuracy of 
DPA is no longer sufficient |4j, and the higher-order corrections in the pole expan- 
sion unlikely can save the situation. Therefore the pinch-technique method and the 
method based on the background-field formalism move forward to foreground, which, 
in principle, can provide the necessary precision (see [5j and [6], and the references 
therein). However, the consecutive application of the mentioned methods implies a 
calculation of a huge volume of additional contributions that formally appear out- 
side the limits of required precision, which is impractical [TJ. So at present hopes 
are pinning on the approach of "complex-mass scheme" (CMS), which avoids men- 
tioned difficulties jU [H]. Nevertheless, in the CMS another problem related to the 
unitarity arises. The point is that the CMS uses the complex-valued renormalized 
masses for unstable particles and this requires an introduction of the complex-valued 
counterterms, which violates unitarity. For this reason the CMS cannot be considered 
as a rigorous procedure [8|. The problem becomes especially topical at calculating 
the contributions in the next-to- next-to-leading order (NNLO). Thus to make the 
calculations up to the NNLO alternative approaches are required. 

A promising candidate for this role is a modified perturbation theory (MPT), first 
proposed in [9] and then elaborated in [10] and For determining the resonant 
contributions the distribution theory is applied in this approach instead of the Dyson 
resummation in whatever form. In essence, the MPT implies a systematic expansion 
in powers of the coupling constant directly of the probability instead of the amplitude. 
This mode allows one to impart the sense of distributions to the propagators squared 
of unstable particles, and on this basis to asymptotically expand the propagators 
squared without the appearance of the divergences in the cross-section. Since the 
object to be expanded (the cross-section) is gauge invariant and the expansion is made 
in powers of the coupling constant, the result of the expansion must automatically 
be gauge- invariant. This implies that the gauge cancellations in the MPT must be 
automatically maintained. Of course, this should be so if the MPT exists, i.e. if it is 
a well-determined method. In the case of pair production of unstable particles this 
property was proved and an algorithm of the calculation of each order of the MPT 
expansion was elaborated [TT] . 

The aim of the given paper is to perform numerical analysis of the convergence 
properties of the MPT series in the case of pair production of unstable particles. 
At once we should notice that in the qualitative sense the outcomes should weakly 
depend on the model under consideration because the choice of a model implies mainly 
a definition of the test function in the presence of which the relevant distributions 
(the propagators squared) are MPT-expanded. So it is reasonable to carry out the 
examination in the framework of a model possessing an exact solution. As such a 
model, we consider the improved Born approximation for the process e + e~ — > 7, Z — > 
ti —> W + b W~b. For simplicity we consider W bosons and b quarks to be stable 
particles, with W being massive and b being massless. At the same time we consider 
realistic corrections to the width of the top quark. This should allow us to get an 
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information about the rapidity of convergence in the realistic case. A similar model 
has actually been considered in [12J by examining the MPT within the next-to- leading 
order (NLO). However, the contribution from the soft massless particles to the process 
of production of unstable particles have been omitted in that work. This was a serious 
omission because the mentioned contributions include Coulomb singularities |13j-[T5] 
appreciably affecting the cross-section. In this paper we improve the calculations of 
[12] (in particular eliminate some bug in the calculations), and carry out numerical 
calculations further up to the NNLO with taking into consideration universal Coulomb 
singular contributions. 

In the next section, we present the basic information about the MPT and detail 
the model in the framework of which we carry out computations. In Sect. [3]we present 
outcomes. In Sect. H] we discuss the results. 



2 MPT and a model for its examination 

The observable cross-section of production and decay of unstable particles, for ex- 
ample in e + e~ annihilation, has the form of a convolution of the hard-scattering 
cross-section with the flux function [TJ, 

a(s)= j ^L<p( s i/s;s)a(s'). (1) 

^min 

Here s is the energy squared in the center-of-mass system, a is the hard-scattering 
cross-sections, <fi is the flux function describing contributions of nonregistered photons 
emitted in the initial state. The s'/s characterizes a fraction of the energy expended 
on the production of unstable particles. For our purposes it is sufficient to take <fi in 
the leading-log approximation. So we put 

4>{z\ s) = (3 e (l - z)^~V - + Z ), e = — (in 4 - lV (2) 

In the case of pair production of unstable particles the double-resonant contri- 
butions are most crucial. Bearing this in mind we write down the hard-scattering 
cross-section in the form 

oo oo 

ct(s) = JJds 1 ds 2 cr(s;s 1 ,s 2 ) (l + 5 c ) . (3) 

s lmin 
■ s 2min 

In this formula <r(s;si,s 2 ) is an exclusive cross-section, 5 C stands for soft massless- 
particles contributions, Si and s 2 are virtualities of unstable particles. In the case 
of the process e + e~ — > 7, Z — > ti — > W + b W~b with massive W and massless b, we 
have Si min = s 2m in = Myy, and s min = 4M^. In a(s ; si, s 2 ) we extract kinematic and 
Breit-Wigner (BW) factors, 

a(s ; si, s 2 ) = — 9(y/s -y/s^-y/s^)^X(s, si, s 2 ) s i> s 2) p{si) p(s 2 ) • (4) 
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Here A(s,si,s 2 ) = [s— (^/~si+^/~S2 ) 2 }[ s — (v^i - V^2 ) 2 ] is the kinematic function, p(si) 
and p(s 2 ) are BW factors. Function $(s; si,s 2 ) is the rest of the amplitude squared. 
Below we consider $ in the Born approximation, and thus only the BW factors are 
subject to the MPT expansion. In the general case, we define the BW factors as 

p(s) = ^ x |A( S )| 2 . (5) 

7T 

Here M is the renormalized mass of the top quark, T is its Born width, A(s) is a 
scalar part of the Dyson-resummed propagator (the spin factor is referred to $), 

A _1 (s) = s-M 2 + ReS(s) + iImS(s), (6) 

ReS(s) and ImE(s) are the real and imaginary parts of the renormalized self-energy. 

In the case of a smooth weight, an isolated BW factor may be represented in 
the form of an asymptotic expansion in the sense of distributions in powers of the 
coupling constant (generating thus the MPT expansion of isolated BW factor). Up 
to and including the NNLO this expansion looks as follows [9]: 

M = S{s -M*) + -Ip^l^—^-^—^j (7, 
+ T, c n(a) { -^-5^(s-M 2 ) + 0(a 3 ). 

Here a is the coupling constant, S(- ■ •) is the 5-function, 8^ n ' is its nth derivative, PV 
means the principal-value prescription. The leading term in (J7]) defines the narrow- 
width approximation. The contributions in the curly brackets appear as a result of 
the naive expansion of the propagator squared; PV makes the poles in this expansion 
integrable. The contributions under the sum-sign correct the contributions of the 
PV poles (in the singular point) so that the expansion becomes asymptotic. The 
coefficients c n (a) are polynomials in a, determined by the self-energy contributions 
of the unstable particle. In an arbitrary UV-renormalization scheme the completely 
explicit expressions for c n (a) may be found in [TO]. In the case of the on-mass- 
shell (OMS) type scheme, they are found in (TTJ. In the latter case the coefficients 
c n within the NNLO are determined by Ii, I2, I3, I[, I", where I n = ImE n (M 2 ), 
/; = ImS^(M 2 ), = ImS^(M 2 ), and by R 2 , R' 2 , where R n = ReS n (M 2 ), R' n = 
ReS^(M 2 ). Here S n is the n-loop self-energy defined in accordance with relation 
£ = a£ 1 + a 2 £ 2 + ---. 

Unfortunately, the weight in our case is not smooth because of the kinematic 
factor in formula (J4]). A solution to this problem is found on the basis of analytic 
regularization via the substitution [A(s, Si, s 2 )] 1 ^ 2 — > [A(s, S\, s 2 )] u - Furthermore, the 
weight $ (1 + 8 C ) may be expanded in powers of Si and s 2 around S\ = M 2 and s 2 = 
M 2 . Then, it becomes possible to analytically calculate singular integrals irrespective 
of details of the definition of the weight. After calculating singular integrals and 
removing the regularization the outcomes remain finite and the expansion remains 
asymptotic [UJ. In principle, this salvages the applicability of the approach, and the 
problem is reduced to numerical calculations only. 



4 



Now we turn to the definition of the model in the framework of which we will 
carry out calculations. At first we notice that within the NNLO the MPT expansion 
of the general BW factor based on propagator §6§ coincides with that based on the 
following "minimal" propagator: 

^nnlo(s) = s-M 2 + aReL l {s) + ia\h + {s-M 2 )I[ + \{s-M 2 ) 2 I'( 



+ a 2 



2 

R 2 + iI 2 + (s-M 2 )R' 2 \ + ia 3 I 3 . 



In fact, after the MPT expansion any contribution not included in (jSJ) appears outside 
the NNLO. However, under the consideration in the conventional-function sense, the 
i?3 and I 2 in the region s — M 2 ~ 0(a) may be assigned to the NNLO, as well. So 
we start from the following modeling propagator: 

a tvW>( s ) = s - M 2 + arReEi(s) + iaImE 1 (s) 



+ a 2 



R 2 + iI 2 + (s- M 2 )(R' 2 + iI' 2 ) +a 3 (R 3 + H 3 ). (9) 



Propagator ([9]) ensures the NNLO precision from the point of view of both the MPT 
and conventional functions. For uniformity we consider ImSi off-shell as well as in 
the case of ReSi. 

Now let us direct our attention to the definition of the two- and three-loop contri- 
butions to propagator (Q. Actually they may be determined by using the only fact 
that they are the on-shell contributions. Specifically, the real parts may be deter- 
mined by basing on the UV-renormalization conditions. But one should remember 
that in the unstable-particles case the OMS scheme may be determined in different 
fashions. In particular, the conventional OMS scheme is determined by the condi- 
tions R n = and R' n = [16J. However, it is inconvenient for the calculations in 
the higher-orders, because the renormalized mass M in this scheme beginning with 
the two-loops is different from the observable mass, and beginning with the two- 
loops generally is gauge- dependent [T7] . The problem is eliminated at considering the 
first renormalization condition in the form M 2 = Re s p , where s p is the pole of the 
propagator, A _1 (s p ) = 0, which means the equating of the renormalized mass to the 
observable mass. The second renormalization condition may be determined by equat- 
ing the imaginary part of the on-shell self-energy to the imaginary part of s p . As a 
result the equality s p = M 2 — i I is established by means of the UV-renormalization 
conditions, where I = ImS(M 2 ). So both the renormalized mass and the imaginary 
part of the on-shell self-energy become gauge- independent. This scheme of the UV 
renormalization was called the OMS scheme in [18] and the "pole scheme" in [19] . In 
this scheme the R 2 , R' 2 and R% are determined as 

R 2 = -hl[ , R' 2 = -hl'l/2 , R 3 = -I 2 I[ - hl' 2 + I 2 R'l/2 . (10) 

The imaginary contributions to the on-shell self-energy are determined, in effect, 
by the unitarity condition. For Ii and I 2 the appropriate relations are 

ah = MT , a 2 I 2 = MaT l . (11) 
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Here To and aT\ are the Born and the one- loop contributions to the width. The /3, 
in the general case, is nontrivially related with the two-loop contribution a 2 T 2 - In 
the OMS scheme this relation is 



a 3 h 



Ma 2 T 2 + rg/(8M) . 



(12) 



Unfortunately, the derivatives of ImE cannot be determined by similar means. 
However the I' 2 , which we need, is a facultative quantity from the point of view of 
the MPT expansion (see above). So we may determine I' 2 by using the approximate 
relationships a 2 ImS2(s) = yfsaT\(s), Ti(s) = Ti x r (s)/r , Tq(s) = almSx (s)/y/s. 
This yields 

-— al x . (13) 

1 n 



a 2 !' 



Now we determine the one-loop self-energy £i(s). In the framework of the model, 
we determine it with contributions of the W boson and b quark only. In this way 
we avoid the IR divergences generally arising at determining ReSi(s). Standard 
calculation in t'Hooft-Feynman gaugqj gives 



a Si(s) 



A(s) 



= A(s 
G F M } 



ReA(M 2 ) - (s - M 2 )ReA'(M 2 



2 

W 



M 2 



B 1 (s]0,M w ) 



1 



(14) 



(15) 



Here -Bi(s;mi,m 2 ) is the Passarino-Veltman function [20] . 

Thus, we have determined all contributions to the propagator ()9]) and thereby 
the BW factors in formula (|3J). Further, by virtue of (j7j) we can determine the MPT 
expansion of the BW factors. The coefficients c n on account of f lTU|) - ffT2"|) . flHJ), ffTB]) 
and HT] follows: 



c 



Cl 



ri 2 

— a — + a 
To 





r 2 " 


r 2 

1 


/r m 2 






8M 2 


\M M 2 



1 w 



M 2 



w , 



0. 



c 2 



-M 2 r 2 



(16) 



Recall that each T n includes an additional factor a, which is conditioned by the vertex 
origin of the width. 

To complete definition of the model, we must determine also the factor (1 + 8 C ) 
in formula ([3]). Let us remember that we have ignored in the self-energy all massless- 
particles contributions that lead to IR divergences. For this reason we have to ignore 
all other soft-massless-particles contributions whose IR-divergent contributions are to 
be cancelled in the cross-section. So, there should remain only the Coulomb singular 
contributions that are not cancelled. Recall that they have the meaning of universal 
corrections arising due to exchanges by soft massless particles (photons, gluons) be- 
tween outgoing massive particles in the limit of small relative velocities. In the case 



1 The gauge independence should be restored after including the higher-order corrections in $, 
and after including the single- and non-resonant contributions in the cross-section 
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of strong-interacting top quarks, it is reasonable to ignore the exchanges by photons 
and to take into account only gluon exchanges. We also restrict our consideration 
to the one-gluon approximation. Then, with taking into account the off-shell and 
finite-width effects, we have [T4J fT5] 



K a q 7T 

Or 



2(3 



\ 2 ( \Pm?-P 

1 arctan 

7T I 2/3lm/3 M 



(17) 



Here k = 4/3 is the group factor, a s is the strong coupling constant, (3 = s 1 y / A(s, s\, s 2 ) 



is the velocity of the unstable particles in the c.m.f., /3m = y 1 — 4(M 2 — iMY)/ s. Fur- 
ther we put T = T in the latter formula. The energy-scale dependence in a s , we take 
into consideration as described in [2Tj . 

So, now the model is completely determined. The cross-section in the model may 
be straightforwardly calculated. We call the result, the "exact" solution. Simulta- 
neously we can calculate the MPT expansion of the cross-section and compare the 
outcome with the "exact" result. Ultimately the expansion should have the form 

cr(s) = ctq(s) + a cti(s) + a 2 a 2 (s) + • • • . (18) 

Here o"o means the cross-section in the LO approximation, a o\ and o?a 2 mean the 
NLO and NNLO corrections, respectively. So, the aoi — &o + oc o\ and 0012 = er + 
a o"i + a 2 a 2 determine the NLO and NNLO approximations. Similarly we denote the 
contributions to the hard-scattering cross-section cr(s). 



3 Results of numerical calculations 

Parameters of the model we determine as follows: M = 175 GeV, M\y = 80.4 GeV, 
and we use the following previously calculated input-data for the width [23J: 

To = 1.56 GeV, 
r + ari = 1.45 GeV, 
r + ari + a 2 r 2 = 1.42 GeV. (19) 

The T = To + aTi + a 2 T 2 we consider as the total width of the top quark. From ffT9~]) . 
we get oTx = —0.11 GeV, and a 2 Y 2 = —0.03 GeV. The a in $(s; si, s 2 ) and <p(z; s), 
we set equal 1/137. All calculations are carried out on the basis of rather general 
FORTRAN code with double precision written in accordance with the formulas and 
instructions described in [TTj . 

In Fig. Hl^a) we present the results of the calculation of the total cross-section a(s) 
above the threshold. The results in percentages with respect to the exact solution 
are shown in Fig. [TJ(b). In the latter figure we place also a result in DPA, where 
<^dpa(s) is determined by the same formulas as in the case of a(s) but by substituting 
$(s;M 2 ,M 2 )(l + 5 c (s;M 2 ,M 2 )) for $(s; s u s 2 )(l + 5 c {s; s u s 2 )) and s-M 2 + iT for 
A _1 (s). The similar results for the hard-scattering cross-section a(s) are presented 
by Fig. |2]^a,b). Let us remember that a(s) is responsible for the distribution over the 
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Figure 1: Total cross-section cr(s). The exact result in the model, we show by thick curve. 
Dotted, short-dashed, and continuous thin curves mean the LO, NLO, and NNLO approx- 
imations in the MPT, respectively. The results are presented in pb (a) and in percents to 
the exact result (b). In panel (b), we show by the long-dashed curve the result in the DPA. 
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Figure 2: Hard-scattering cross-section <r(s). The notation is the same that in Fig. [TJ 



invariant mass of the tt system and therefore is of interest, as well [21]. In Fig. [3] we 
show the results separately for the NLO and NNLO corrections to a. In Tabled] the 
results are represented in the numerical form at the characteristic energies accessible 
at the planned e + e~ colliders. In the last column the numbers in parenthesis represent 
the uncertainties in the last digits. (See discussion of their determination in |25j.) In 
the other columns the uncertainties are omitted as they appear in the digits that are 
not shown. In the lower positions in the Table the results are presented in percentages 
with respect to the exact result in the model. 

The above outcomes exhibit very stable behavior of the NLO and NNLO approx- 
imations in the energy region beginning with approximately 400 GeV. (In this region 
simultaneously the right hierarchy of the corrections is established, <jq < o\ < 
The accuracy of the NNLO approximation is established greatly high in this region. 
In particular, at 400 GeV < y/s < 600 GeV it is within ±0.5%. At increasing energy 
the accuracy in relative units is slightly decreasing, but the cross-section is decreas- 
ing, too, so that the effective precision of the description remains approximately the 
same (because the ratio of the discrepancy to the quantity y/a, which characterizes 
statistical error, is approximately constant). In contrast to the above picture, the 
DPA exhibits greatly unstable behavior; its discrepancy varies from +3.0% to -7.4% 
in the energy region 400 GeV < \/s < 1500 GeV. On the whole, such a behavior 
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Figure 3: Corrections o\ and o 2 (dashed and continuous curves, respectively). 



Table 1: The results of the calculation of the total cross-section in pb. 
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0.04 



0.02 



coincides with the expected one for DPA, including the order of magnitude of the 
discrepancy in the Born approximation for <3> [2J. 

To conclude this section three important remarks are in order. First we note 
that the results expressed in relative units are almost insensitive to the choice of the 
test function $. In particular, the turning-on/off of the Coulomb factor has very 
small effect. For instance, at yfs = 500 GeV this leads to 0.7%-modification of the 
ratio 0012/0" and at y/s = 1500 GeV does less than 0.1%. (Although, the variation 
of the absolute value of the cross-section is considerable in both cases: about 35% 
and 20%, respectively.) The second remark concerns a large value of the correction 
cr 2 in comparison with the discrepancy — cr 012 . For example, at y/s = 500 GeV 
they constitute 5.3% and 0.4% of 0, respectively. However we think that this is an 
incidental unbalance as 02 gains its value mainly due to the correction a 2 T 2 to the 
width, which in our case exhausts the corrections, and simultaneously a 2 r 2 is quite 
large (approximately 2% of T). If we put everywhere a 2 T 2 = 0, then at y/s = 500 GeV 
the o"2 decreases to 2.1% with the discrepancy remaining within the 0.5%-interval. On 
the other hand, if we put ct 2 r2 = only at calculating the coefficients c n without the 
change of the model itself, then the 2 becomes almost the same as in the latter case, 
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but the discrepancy increases to 3.6%. The third remark concerns the ill-convergent 
property of the MPT in the near-threshold region. With the energy approaching the 
threshold the accuracy and stability of MPT rapidly become worse. This manifests 
itself in the violation of the hierarchy of the corrections and then in the blowing up 
of the corrections. Actually, this behavior was predicted in [TTJ. To prevent this 
difficulty another mode of the MPT near threshold must be applied [TTJ that implies 
Taylor expansion of a(s) both in powers of a and in powers of s — 4M 2 , where 4M 2 
is the threshold. An alternative method implies a secondary Dyson resummation in 
the framework of the MPT approach [TUl 122] . 

4 Discussion 

Although above calculations have been carried out in the framework of a model, the 
obtained outcomes expressed in relative units to a large extent are model- independent. 
By this we mean that the outcomes are weakly sensitive to the choice of the test 
function determined by the model under consideration. We verified this property by 
carrying out calculations with various test functions and found that the influence of 
the test function manifests itself mainly in a factor common for different contributions 
to the cross-section. In particular, even very large variation in the test function that 
appear at the turning-on/off the Coulomb factor, in relative units leads to small 
modifications of the outcomes. 

On this basis we suppose that the loop corrections to the test function will lead 
in relative units to small modifications of the outcomes, too. In particular, our re- 
sult about the 0.5%-accuracy of the NNLO approximation near the maximum of the 
cross-section, should remain in force at turning-on the loop corrections. Moreover, 
one can further improve the results if applying the MPT on the background of the 
loop corrections only, and considering the Born contribution in the old fashion with 
the Dyson-resummation in the unstable-particles propagators — on analogy of ac- 
tual practice of application of DPA [HE]- In this case the discrepancy in the MPT 
description will be diminished by a factor 0(a). 

Another aspect of the problem of model-dependence of our results concerns the 
corrections to the width of unstable particles. Recall that these corrections determine 
coefficients c n , which are crucial for the definition of MPT expansion. We have 
considered the case with rather large corrections to the width (7% and 2% in the NLO 
and NNLO, respectively). At diminishing these corrections, the MPT corrections to 
the cross-section should diminish, too. At least, we have observed this property in 
the framework of the model under consideration. On this ground we can expect the 
improving of the precision of description at transiting from the top quarks to EW- 
only interacting particles, for instance to the VT-bosons, because the corrections to 
the width are lesser in the latter case. 

As regards the application of our results to the description of realistic processes 
with the top-quark pair production, we should stress that our calculations simulate the 
main contribution to the cross-section as they cover the double-resonant contributions. 
So on the basis of our results we can judge about the precision that must be achieved in 



10 



realistic calculations. Fortunately, the accuracy of the NNLO approximation detected 
in our analysis, is satisfactory from the point of view of the ILC requirements. Really, 
assuming that at the ILC several hundred thousands of the tt events is expected, we 
conclude that the calculation of the cross-section is needed with a few per mille 
accuracy. As we have seen above, this, in general, is ensured by the NNLO in the 
MPT. 

In summary, we have shown that the MPT stably works at the energies near the 
maximum of the cross-section and above at the description of the total cross-section 
for the pair production and decay of fundamental unstable particles. We have found 
also that in the mentioned energy region the MPT provides very good precision 
within the NNLO. In particular, at the ILC energies in the case of the top-quark 
pair production the NNLO approximation provides 0.5%-precision of the description. 
The further increase of the precision is possible at the proceeding to the NNNLO, 
possible on the basis of the results of [TT], or at the proceeding to the compound use 
of the MPT, when the loop corrections are treated completely in the framework of 
the MPT while the Born contribution to the cross-section is taken into consideration 
in the old fashion with the Dyson resummation in the unstable-particles propagators. 
On the whole, the MPT method is a real candidate for carrying out high-precision 
calculations needed for ILC. 
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